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ABSTRACT 

The reconstruction of the peculiar velocity field from the 1.936 Jy /_R/l 5 selected sample 
of galaxies is compared to a similar reconstruction from an optically selected sample. A 
general method for combining different samples to reconstruct a self-consistent density 
and peculiar velocity field is presented. The method is applied to determine how sensitive 
the derived peculiar velocity field is to the characteristics of the sample used. The 
possibility that the IRAS galaxies do not trace the general galaxy population is explored 
adopting a simple model of linear biasing between the IRAS and optical samples. We 
find that the velocity fields derived from the two samples are consistent, within the 
estimated shot noise error, for the case of no relative bias. This result suggests that the 
predicted peculiar velocity field based on /_R/l 5 samples is not sensitive to the sampling 
properties of IRAS galaxies. Combined with previous suggestion of a relative biasing 
of IRAS galaxies on small scales (about bh~^Mpc), this result suggests scale dependent 
biasing. 

Key words: large scale structure of the universe - cosmology: observations - galaxies: 
clustering - galaxies: distances and redshifts 



1 INTRODUCTION 

Under the assumption that galaxies trace the matter distri- 
bution, it is possible to reconstruct a self-consistent density 
and peculiar velocity field from observed redshifts. Compar- 



ison of the predicted peculiar velocity field and the avail- 
able data from various redshift-distance surveys is a power- 
ful approach to the study of large-scale structure, which can 
be used to test the gravitational instability picture for the 
growth of structures, to measure the relative distribution of 
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luminous and dark matter and to determine the quantity 
Q /h, where Q is the density parameter and h is the bias 
of the galaxy distribution relative to the underlying matter 
density (e.g. Peebles 1980, Dekel et al. 1993 and references 
therein) . 

A basic requirement for the self-consistent reconstruc- 
tion procedure is the availability of a uniform "whole-sky" 
sample with complete redshift information. This require- 
ment has led to the extensive use of samples drawn from 
the IRAS catalog of galaxies, taking advantage of its wide- 
angle coverage and uniformity. In particular, this method 
has been applied by several authors (Yahil 1988, Strauss & 
Davis 1991, Strauss et al. 1992) to a sample of galaxies with 
complete redshift information to a limiting flux of 1.936 Jy 
at 60jUm in the IRASVAC database. This sample covers over 
87% of the sky. 

Unfortunately, the IRAS saTa^Ae is rather dilute, delin- 
eating poorly some of the major structures of the galaxy 
distribution as compared to those observed in optical sur- 
veys. In addition, it is well known that the IRAS galaxies 
represent a special class of galaxies and are relatively de- 
ficient in high-density regions, as most IRAS galaxies are 
late-type spirals. Some authors have suggested that the dis- 
crepancy between the optical and IRAS dipoles computed 
from the two-dimensional distribution could imply that the 
distribution of //? A 5 galaxies is biased relative to the optical 
galaxies (e.g. Lahav, Rowan- Robinson & Lynden-Bell 1988). 

Here we examine how sensitive the predicted flow is to 
these effects. This is done by comparing the velocity field ob- 
tained using the //?A 5 sample with that derived from a com- 
bined Optical-IRAS saTa^Ae. The combined sample used here 



takes advantage of the denser sampling offered by optical 
redshift surveys over roughly 50% of the sky, and the wide- 
angle coverage of the IRAS survey. We use this combined 
sample to investigate the degree of relative bias between the 
distribution of optical and IRAS galaxies. We also use the 
larger optical sample to estimate the error in the predicted 
velocity field due to shot noise. The possibility of a relative 
bias between optical and IRAS galaxies has recently been 
investigated by Strauss et al (1992) by direct comparison of 
the smoothed density fields obtained from the IRAS sara^Ae 
and the CfAl, SSRS and SPS optical samples (Huchra et 
al 1983, da Costa et al 1991, Dressier 1991). They show 
that on large scales the two density fields are consistent 
with them being drawn from the same parent distribution. 
However, on smaller scales (k, 500 km s~ ) there are sta- 
tistically significant differences between optical and IRAS 
galaxies. Since a range of scales contribute to the peculiar 
acceleration of a galaxy it is unclear how this may affect 
the predicted peculiar velocity which is frequently used for 
comparison with those measured using redshift independent 
distance estimates. Here we investigate what effect these dif- 
ferences may have on these predictions. 

Over the years several attempts have been made to cal- 
culate the Local Group velocity based on optical samples 
using different approaches to treat the low galactic latitude 
zone (e.g., Davis & Huchra 1982, Pellegrini & da Costa 1990, 
Hudson 1993). There are several important differences be- 
tween the current work and previous approaches, which in- 
clude the optical sample used and a more realistic interpo- 
lation across the obscuration zone by utilizing the informa- 
tion available from the //?A 5 catalog. More importantly, this 



analysis is not restricted to the Local Group but considers 
the global characteristics of the velocity field. 

In section 2 we discuss the details of the method we 
have utilized for computing the peculiar velocity field, the 
construction of the combined Optical- IRAS sample and how 
errors were estimated. In section 3 the results of different 
models are presented together with a discussion of their im- 
plications. A brief summary of our main conclusions is pre- 
sented in section 4. 



2 THE METHOD 



2.1 Physical Basis 



The measured radial velocity of a galaxy is the sum of its 
real distance and peculiar velocity 



V = r -\- u ■ r , 



(1) 



where v is the measured radial velocity in the Local Group 
rest frame, r is the distance from the sun, u is the peculiar 
velocity vector, and r is the unit vector pointing from the 
observer along the line of sight. Since only relative distances 
are relevant to this discussion, we have expressed distances 
in units of km s~ in order to avoid specifying a value for 
the Hubble constant. 

If the peculiar velocities are gravitationally induced and 
galaxies trace matter, they will be related to fluctuations 
of the galaxy distribution. Linear theory of the growth of 
density perturbation predicts that the peculiar velocity of 
any given galaxy is 
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where g is the peculiar acceleration of the galaxy. It is calcu- 
lated from the potential (f>, obtained by numerically solving 
Poisson's equation 
3. 



V",^= -Q(1-F5„) 



(3) 



2g 

3Q"4 



(2) 



where &rn = Pmj < Pm > — 1 is the matter density fluctua- 
tion; <> indicates an average over the sampled volume. 

Once the distribution of galaxies is specified in redshift 
space, the peculiar acceleration g and the peculiar velocity 
u can be computed from equations (3) and (2). The pecu- 
liar velocity can then be used in equation (1) to improve 
the estimates of the positions for the galaxies with mea- 
sured redshift v, iteratively. The method breaks down in 
regions of high density, where there is no unique distance 
for a given redshift. In these regions, we follow method 2 of 
Yahil et al (1991) and adopt the distance which is the aver- 
age of the minimum distance where the computed redshift 
is 200 km s~ less than the observed redshift and the max- 
imum distance where the computed redshift is 200 km s~ 
more than the observed redshift. For every galaxy, the pe- 
culiar velocity was computed at 100 distances along the line 
of sight and for each such point the corresponding redshift. 
From the resulting redshift-distance relation, galaxies which 
need this special treatment were identified and their distance 
computed. A different strategy to treat such regions has been 
discussed by Freudling et al. (1991). 

Since the computation of the peculiar velocities relies on 
a linear approximation, it cannot be applied within clusters. 
Therefore, clusters from the list given by Yahil et al. (1991) 
have been collapsed in redshift space by setting the initial 
distances and positions of the cluster galaxies to that of the 
center of the cluster. 
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While the physical basis of the algorithm used here to 
compute the accelerations g is identical to the one used 
by Yahil et al. (1991), the numerical approach is different 
(Freudling et al. 1991). Instead of adding the contribution of 
individual galaxies, the gravitational acceleration g is com- 
puted from the smoothed density field with a particle-mesh 
code via equation (3). The density field is calculated on a 
grid with a 500 km s~ cell size. The total dimension of the 
grid is 64 , i.e. the sphere with radius 8000km s~ on which 
& is not equal to zero occupies only a fraction of the cube. 
These dimensions were chosen so that galaxies do not inter- 
act with the mirror image created by the periodic boundary 
conditions of the particle-mesh code (see Freudling et al., 
1991). The weight of each galaxy is assigned to the grid 
so that the resulting smooth distribution does not intro- 
duce any dipole or quadrupole moments (TSC assignment 
scheme). The weight assigned to galaxies is the inverse of 
the appropriate selection function. Davis & Huchra (1982) 
have shown that this assignment is virtually identical to the 
minimum variance estimator. Finally, the relative density 
fluctuation & is computed for grid points in the surveyed re- 
gion (v < 8000 km s~ ). The advantage of this algorithm is 
that it provides direct access to the & used as input to the 
linear perturbation theory. We take advantage of this for the 
interpolation over the galactic plane and for the implemen- 
tation of the bias of the //?A 5 sample relative to the optical 
sample. In principle, a simple nonlinear correction, as sug- 
gested by Nusser et al. (1992), could also be easily built into 
the model. 



2.2 Construction of the Density Field 

2.2.1 Samples 

In this paper we have considered two samples: the first is the 
IRAS 1.936 Jy sample of Strauss et al. (1992); the second 
one is a combination of the optical sample available in the 
galactic caps and the part of the IRAS sample located at 
low galactic latitudes, which are not covered by the former. 
The optical sample consists of galaxies brighter than 
"^B(o) =14.5. In the northern galactic cap (b>40°) we use the 
CfAl sample (Huchra et al. 1983). In the southern galactic 
cap (b < -30°) we use the so-called SGC sample assembled 
by Pellegrini et al. (1990) from different catalogs (CGCG, 
FSO, MCG) with redshifts taken primarily from the SSRS 
(da Costa et al. 1991) and the equatorial survey (Huchra 
et al. 1993). This combined catalog has been shown to be 
fairly homogeneous and its magnitude scale to be statis- 
tically consistent with the B(0)-Zwicky system adopted in 
the CfAl survey. Recently, improvements have been made 
to this catalog south of 5 = —17.5° and below h = —30° by 
adopting the Bt magnitudes listed in the FSO photomet- 
ric catalog (Lauberts & Valentijn 1989), converted to the 
B(0) system, instead of using magnitudes estimated from 
diameter-magnitude relations. A more detailed discussion 
on the accuracy of the resulting magnitudes has been pre- 
sented by Alonso et al. (1993). 



2.2.2 Density Field for the Optical-IRAS Sample 

If only one sample is considered the density fluctuation S 
can be computed directly from the number counts per vol- 
ume p as S = p/ < p > —1. This implicitly assumes that the 



surveyed volume represents a fair sample of the universe and 
that the average density fluctuation < 5 > is zero within this 
volume. In order to impose the same assumption that the 
total volume averaged S of the combined Optical- IRAS sam- 
ple equals zero, the differences in the mean density between 
the two distinct samples must be taken into account. 

This can be achieved in the following way. To compute 
S for our combined sample from the number densities po and 
Ps in the optical and IRAS part of the sample, respectively, 
we use the substantial overlap between these samples in the 

galactic caps (h < —30° or 6 > 40°) to define 
' Po in caps 



(4) 



Pt 



<Po>o 
<P,>a 



in plane 



where <>c is the average of a quantity in the galactic caps. 

With this expression, the density p is computed over the 

entire volume, and its average used to determine the density 

fluctuation 

' Po/ < p > —1 in caps 



(5) 



p. <Po>c 



1 in plane 



<P> <P,>a 

In order to account for the possibility of a relative bias 
between the IRAS and optical samples we introduce the rel- 
ative biasing factor 6r = &optl&iras, where &,ras and &opt are 
the density fluctuations of the IRAS and optical samples, 
respectively. Here we have adopted a simple linear biasing 
scheme between these samples, assuming that 6r is a con- 
stant independent of position. For models which utilize the 

combined sample, equation (5) changes to 
' Pol < p' > —1 in caps 

(6) 
where the density < p > is adjusted to make < 5 >= as 
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before. In the next section we use this equation to compute 
S for different values of 6r . 

For grid points at high redshifts, beyond the velocity 
limit imposed on our samples, the density fluctuation S was 
set to zero. At low galactic latitudes, we follow Yahil et al. 
(f99f) and interpolate over the region |6| < 5° by setting 
the density equal to the average of the density in the 2 
strips 5° < |6| < f5° in fOOO km s"^ redshift bins and fO° 
longitude bins. This interpolation can again contribute to 
the mass contained within the sphere defined by the maxi- 
mum redshift, because the average density fluctuation S in 
the interpolation zone is not necessarily zero. This contribu- 
tion has to be taken into account when the average density 
< p > is computed. Once the initial density field has been 
constructed the system is allowed to evolve as outlined be- 
low. 



2.3 Implementation 

The complete computational procedure was as follows. As 
a preparation for the iterative part, the observed redshifts 
of the galaxies were corrected to that relative to the center 
of the Local Group, which was assumed to move toward 
I = 270°, 6 = 0° with 300 km s'K The initial distance of 
all field galaxies was set to their corrected redshifts, and 
clusters collapsed as described above. 

Fach iteration started with the computation of the se- 
lection function, following the procedure adopted by Davis 
& Huchra (f982). For the combined Optical-IRAS samples, 
the function for the IRAS part of the sample was assumed to 
be that derived from the corresponding all-sky IRAS model 
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in order to reduce the uncertainty due to the smaller sky 
coverage. 

To improve the convergence of the procedure, the den- 
sity parameter Q was set to a fraction of its final value in 
the first iteration, and gradually increased to its final value, 
always taken to be unity. This has the effect that the den- 
sity field, and therefore the gravitational field, changes only 
slightly between iterations. 

The next step is to construct the density field from the 
discrete distribution of the galaxies. This includes the distri- 
bution of the galaxy weights to the grid, biasing and inter- 
polation over the galactic planes which has been discussed 
in detail in previous section. Subsequently, equation (3) is 
solved to compute the gravitational field and peculiar veloc- 
ities by equating the matter density fluctuation &rn to the 
optical galaxy density fluctuation given by equation (6). To- 
gether with our assumption that Q = 1 this is equivalent to 
models with /3 = Q /6 = 1, where h is the bias factor of 
the optical sample relative to the matter distribution. 

The distances of the galaxies is adjusted for the next 
iteration to be the average between their previous distance 
and the distance which produces the observed redshift given 
the computed peculiar velocity, or using the procedure de- 
scribed above where necessary to avoid ambiguities. 

The rms of this distance correction as well as the 
maximum adjustment was recorded, and the next iteration 
started with the new distances. In typical runs, the models 
converged to an rms of about 10 km s~ in 5 iterations after 
the final set of parameters was reached. 

A preliminary value for the relative weight between the 
IRAS and the optical part of the sample < po > / < Pi > 



was computed from the real space distributions for the IRAS 
and optical models restricted to the galactic caps (see dis- 
cussion below) and was found to be 2.40. Subsequently, this 
value was used as input for a preliminary version of the 
whole-sky, unbiased IRAS and Optical-IRAS models, and 
was recomputed from the real space distribution of these 
models, yielding < Po > I < Pi >= 2.45. This value was 
adopted in all other models. This factor computed sepa- 
rately for the northern and southern cap resulted in values 
which differ from each other by 22%. 

2.4 Error Estimates 

In order to estimate the shot noise typical of IRAS-Yike 
samples we have used the larger Optical-IRAS sample in 
the following manner. Since most of the contribution to the 
shot noise comes from large distances, we created subsets of 
the large optical sample which had approximately the same 
number density as the IRAS sample near our velocity limit 
of 8000 km s~ . The optical sample at 8000 km s~ contains 
about twice as many galaxies per volume as the IRAS sara- 
ple. Therefore we randomly divided the optical sample into 
two. The purpose of creating these diluted samples was to 
examine the difference between the IRAS and the Optical- 
IRAS velocity models simply due to the sampling, i.e. the 
shot noise. These diluted samples contain about twice the 
number of galaxies per solid angle as the IRAS sample in 
the galactic caps, which reflects the different radial distribu- 
tion of the two samples. These optical subsamples were then 
combined with the IRAS sample restricted to the galactic 
plane in the same way as it was done for the complete optical 
sample. 



The difference in the predicted pecuUar velocities of the 
two diluted samples was used as an empirical estimate of the 
shot noise as a function of distance. Local Group velocities 
derived from the models without biasing are given in Ta- 
ble 1, where the number of IRAS and optical galaxies in the 
diluted samples and the direction and amplitude of the Lo- 
cal Group velocity vector predicted from these models are 
given. Figure 1 quantifies the difference between those mod- 
els as a function of distance. The upper panel shows the 
average difference of the absolute value of the peculiar ve- 
locities as a function of distance. This comparison includes 
both galaxies within and outside of the tripled valued re- 
gions. Galaxies have been binned in 500 km s~ distance 
bins. The error bars show the rms / y N in each bin, where 
A^ is the number of galaxies in that bin. The middle panel 
shows the angular difference between the velocity vectors 
in both models. The binning and error bars are computed 
as in the upper panel. The lower panel shows the distribu- 
tion of the differences in the x, y and z components of the 
peculiar velocities. Note that this distribution is not nec- 
essarily expected to be symmetrical because differences are 
computed at the positions of galaxies, which might sample 
the velocity field preferentially where differences are nega- 
tive or positive. 

Errors in the amplitude are less than 200 km s~ , and 
angular differences are in the range of 20° to 30°. As can be 
seen in figure 1, the shot noise is small at small distances. 
This is because velocities are relative to the Local Group 
rest frame and therefore by definition zero at the origin in 
all models, and the sampling is densest at small distances. 
The shot noise increases with distance as the sampling be- 
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comes more dilute. At distances greater than the velocity 
limit of the sample (8000 km s~ ), all contributions to the 
gravitational acceleration come from the inner volume and 
therefore the shot noise decreases as this limit is approached. 
Since all peculiar velocities are computed in the Local Group 
rest frame, predicted peculiar velocities in each model con- 
verge towards the projection of the predicted Local Group 
velocity for galaxies at large distances. Therefore, for very 
large distances, the differences between any two models ap- 
proach asymptotically the difference in the estimated Local 
Group velocity for each sample. 

This comparison was used as an estimate of expected 
differences due to shot noise between the full Optical-IRAS 
model and the IRAS model. However, this error estimate 
was based on two IRAS-Yike diluted samples. Due to the 
denser sampling of the optical distribution, this will be an 
overestimate. To take this into account, we divide the orig- 
inal error estimate by v2, which will be used as a lower 
bound to our error estimate, whereas the original error will 
be used as an upper bound. 

3 RESULTS 

First, to test our procedure, the velocity field was com- 
puted using only the IRAS sample, limited to a redshift 
of 8000 km s~ and taking 6r = 1. For this case the pre- 
dicted velocity field, in cartesian supergalactic coordinates, 
is shown in figure 2. In this representation the SGX-axis 
points to SGL=0° and the SGY-axis points nearly in the 
direction of the north galactic pole. Fach row contains three 
parallel and non- overlapping 3000 km s~ thick slabs, with 
the middle slab centered at the origin. In these figures the 
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SGY=0 plane is almost coincident with the galactic plane 
while the SGZ=0 plane is coincident with the supergalactic 
plane. The unit throughout is km s~^ and the amplitude 
of the vectors are in the same units as the grid. The data 
is presented in the Local Group rest frame. In each panel 
the vectors represent the projection onto the corresponding 
plane of the predicted three-dimensional velocity vector. The 
velocity field of figure 2 shows good agreement with those 
of Yahil (1988), who plotted galaxies within a 45° wedge 
about the Supergalactic plane, confirming that the results 
are insensitive to the details of the method used. This can be 
seen, for instance, in the special case of the induced velocity 
of the Local Group for 6r = 1, given in Table 2. The columns 
are as follows: column (1) the sample used; column (2) the 
bias factor; columns (3) and (4) the number of IRAS and 
optical galaxies in the sample; column (5) the number of 
galaxies in triple valued regions in the model, columns (6) 
and (7) the galactic longitude and latitude of the direction 
of the Local Group velocity vector; column (8) the ampli- 
tude of the Local Group velocity and column (9) the angle 
between the computed Local Group velocity vector to the 
cosmic microwave background dipole (Smoot et al. 1992). 
The predicted Local Group velocity for the unbiased model 
is within 4 degrees in direction and 5% in amplitude of that 
obtained by Strauss et al. (1992) for a similar model. 

Next, several models using the combined Optical-IRAS 
samples were computed, for values of relative biasing fac- 
tor 6r in the range 0.7-1.8. For comparison, similar models 
with the IRAS saTa^Ae were computed. These IRAS raoAeXs, 
with different values of 6r, are intended to simulate whole- 
sky optical samples within the framework of a linear biasing 



model. The velocity field derived for 6r = 1 is illustrated 
in figure 3. The velocity vectors are only shown at the po- 
sition of the IRAS galaxies. The similarities between the 
IRAS and Optical-IRAS samples are striking. In order to 
assess their differences, the velocity fields derived from the 
Optical-IRAS sara^Ae were subtracted from those derived by 
the corresponding biased IRAS models. The residual fields 
are shown and figures 4 and 5. Comparison of these figures 
clearly illustrates that the difference between 6r = 1 IRAS 
and Optical-IRAS models are much smaller than between 
the 6r = 1.3 models. In figures 6 and 7 the rms differences 
as a function of distance were compared with the estimated 
errors obtained from the diluted samples. As can be seen 
from the examination of these figures for the case 6r ~ 1 the 
departures are small while for 6r = 1.3 significant departures 
can be observed out to 4000 km s~ . 

Instead of looking at the global velocity field we can 
focus our attention on the peculiar velocity induced to the 
Local Group. In Table 2 we list the values obtained for the 
Local Group velocity for some of the different models de- 
scribed above, while in figure 8 we show the vector difference 
of the Local Group velocity predicted by the Optical-IRAS 
and IRAS sample as a function of the relative biasing fac- 
tor. It is interesting to note that the Local Group velocity 
vector for the Optical-IRAS and IRAS models with 6r ~ 1 
is essentially the same. This result could indicate that no 
linear relative bias between optical and IRAS is required. 
Another possibility is that the velocity field even at large 
distances from the galactic plane is primarily determined by 
structures near it, where our combined Optical-IRAS sam- 
ple is by construction identical to the IRAS sample. This 



explanation would be surprising since major structures such 
as the supergalactic plane, Virgo and the Southern Wall are 
located at high galactic latitudes. We tested this hypothesis 
by running the same models with the density fluctuations S 
in the galactic plane set to zero. The Optical-IRAS models 
are thus reduced to the optical sample only and the IRAS 
sample is restricted to the same volume. The resulting Lo- 
cal Group velocities for 6r = 1 are also given in Table 2. 
The agreement between the Local Group velocity vector de- 
rived from these two models is again excellent both in direc- 
tion and amplitude, although they are now constructed from 
completely different galaxy samples. We conclude that any 
linear relative biasing between optical and //? A 5 galaxies is 
small. 
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we are directly comparing independent optical and IRAS 
samples. These results suggest that the predicted velocity 
field derived from the //?A 5 catalog is robust and that IRAS 
galaxies are fair tracers of the galaxy population. The good 
agreement of the velocity fields from these two samples also 
confirm, a posteriori, our original claim that it is possible to 
construct a fairly homogeneous "all-sky" optical sample. 
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4 CONCLUSIONS 

In this paper we have compared the peculiar velocity field 
predicted from an IRAS and an optical sample in order to 
directly test the effects that diluteness and undersampling 
of high-density regions may have on the computed velocity 
field from the IRAS catalog of galaxies. Since the available 
optical samples have limited sky coverage, we have adopted a 
combined Optical-IRAS catalog and modeled the differences 
between the density fluctuations in terms of a linear biasing 
model. 

We find that the predicted flow is similar, within the 
estimated shot noise, if the relative bias is small. This result 
suggests that the relative biasing between the IRAS sample 
and the optical sample is small on scales which induce the 
peculiar velocities. We find that this is true even when the 
samples are restricted to the galactic caps, in which case 
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Table 1. Diluted Models. 

n n (kms-i) 

970 1945 241 38 734 

967 1937 246 32 749 

Table 2. Predicted Local Group Velocities. 



sample K n^ras rtopt "trtp hg big 



Vlg 



n n (kms-i) (°) 



IRAS 0.7 1897 

IRAS 1.0 1893 

IRAS 1.3 1884 



7 251 36 

235 252 34 
620 253 35 



IRAS in caps 1.0 928 

optical in caps 1.0 3851 



574 23 

726 21 

881 21 



Optical-IRAS 0.7 970 3873 99 241 39 646 31 

Opttcal-IRAS 1.0 969 3873 244 242 34 719 29 

Opttcal-IRAS 1.3 961 3932 499 249 35 776 24 



232 65 509 44 

207 59 501 54 
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Figure Captions 

Figure 1- Differences in the velocities fields of the two di- 
luted models. The upper panel shows the differences 
in the absolute value of the predicted velocities as a 
function of distance. Galaxies have been binned into 
500 km s~^ bins according to the average distance in 
the two models. The error bars shown are the rms di- 
vided by the square root of the number of galaxies in 
each bin. The middle panel shows the differences in 
the direction of the velocity vectors in degrees, again 
binned into 500 km s~ bins. The error bars are com- 
puted as in the previous panel. In the lower panel the 
distribution of the differences in x, y and z direction 
are shown. 

Figure 2- Predicted velocity field computed from the den- 
sity field of IRAS galaxies. Fach panel shows the veloc- 
ity components, projected onto the indicated planes, 
of galaxies within 3000 km s~^ slabs. To avoid over- 
crowding, velocity vectors are shown only for a sample 
semi-volume limited at 4000 km s~ . Distances and ve- 
locities are both in km s~^ . Panels are ordered from 
the top-left (a) to the bottom-right (i). We use super- 
galactic coordinates as explained in the text. 



computed from the density field of the Optical-IRAS 
and IRAS sa.TaY>\e. The relative biasing 6r was 1 in both 
models. 

Figure 5- Same as figure 4, for 6r = 1.3. 

Figure 6- Differences in the velocities fields between the 
Optical-IRAS and IRAS models for the case 6r = 1. 
The meaning of the three panels are described in the 
caption to figure 1. The solid and dashed lines repre- 
sent the upper and lower bounds of our error estimate. 

Figure 7- Same as figure 6 for the case 6r = 1.3. 

Figure 8- Differences in the Local Group velocity vectors as a 
function of the relative biasing factor. The upper panel 
shows the fractional change of the amplitude, while the 
lower panel shows the difference in the angular direc- 
tion. 



Figure 3- Same as figure 2, showing the predicted velocity 
field, at the position of the IRAS galaxies computed 
from the density field of the Optical-IRAS saTa^Ae and 
6r = 1. 



Figure 4- Same as figure 2, but showing the difference of the 
predicted velocity, at the position of IRAS galaxies. 



